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TRPM6 is crucial for human Mg 2+ homeostasis as patients carrying TRPM6 mutations develop 
hypomagnesemia and secondary hypocalcemia (HSH). However, the activation mechanism of TRPM6 has 
remained unknown. Here we demonstrate that phosphatidylinositol-4,5-bisphophate (PIP2) controls 
TRPM6 activation and Mg 2+ influx. Stimulation of PLC-coupled Ml -receptors to deplete PIP 2 potently 
inactivates TRPM6. Translocation of over-expressed 5-phosphatase to cell membrane to specifically 
hydrolyze PIP 2 also completely inhibits TRPM6. Moreover, depolarization-induced-activation of the 
voltage-sensitive-phosphatase (Ci-VSP) simultaneously depletes PIP 2 and inhibits TRPM6. PLC-activation 
induced PIP 2 -depletion not only inhibits TRPM6, but also abolishes TRPM6-mediated Mg 2+ influx. 
Furthermore, neutralization of basic residues in the TRP domain leads to nonfunctional or dysfunctional 
mutants with reduced activity by PIP 2 , suggesting that they are likely to participate in interactions with PIP 2 . 
Our data indicate that PIP 2 is required for TRPM6 channel function; hydrolysis of PIP 2 by PLC-coupled 
hormones/agonists may constitute an important pathway for TRPM6 gating, and perhaps Mg 2+ 
homeostasis. 



Magnesium (Mg 2+ ) is the most abundant divalent cation in the mammalian cell and is essential for 
numerous fundamental cellular processes, including cell cycle, channel regulation, ATPase activity, 
metabolic regulation, and various signaling pathways 1 . Mg 2+ deficiency has been implicated in many 
diseases, ranging from neurological to cardiovascular diseases 2, 3 . Mg 2+ homeostasis is therefore tightly controlled 
by maintaining the equilibrium between intestinal Mg 2+ absorption and renal Mg 2+ excretion/re- absorption. 

Several Mg 2 * transporters and channels have been implicated to be important in Mg 2+ absorption and/or 
reabsorption 411 . Most notable, loss of function of TRPM6 causes familial hypomagnesemia and secondary hypo- 
calcemia (HSH) 8, 10 . TRPM7, the closest homologue of TRPM6, was demonstrated to be essential for cellular Mg* + 
homeostasis in DT-40 lymphocytes 11 . However, it's role in Mg 2+ homeostasis in mice remains disputed 12, 13 . TRPM6 
and TRPM7 are unique bi-functional channel proteins with protein kinase function 1418 . TRPM7 is ubiquitously 
expressed in various cells and tissues and has well defined functions in embryonic development 12, 13, 19 , neuronal cell 
death 20 , and a variety of other functions 2124 . Different from TRPM7, TRPM6 expression is restrained in the 
epithelial cells in intestine and distal convoluted tubule (DCT) of nephron 9, lu , consistent with its central role in 
controlling Mg 2+ homeostasis. 

Mg 2+ homeostasis is regulated by a variety of hormonal and pathological conditions 1 . As a gatekeeper of 
human Mg 2+ homeostasis 25 , TRPM6 has been demonstrated to be regulated at expression levels 2, 26 by hormones 
such as estrogen 26 and Angll 27 , metabolic acidosis/alkalosis 28 , immunosuppressant tacrolimus 29 , diuretics 
Thiazide 30 , and EGF 31 . However, the gating mechanism of TRPM6, the key property which controls Mg 2+ influx, 
has remained elusive. Like TRPM7, TRPM6 is inhibited by millimolar concentration of intracellular Mg 2+ 
([Mg 2+ ]j); therefore, it only constitutively opens to a small degree under physiological [Mg 2+ ]i 18 . Both TRPM6 
and TRPM7 are permeable to Ca 2+ and Mg 2+ under physiological pH, and conduct monovalent Na + currents at 
acidic extracellular pH 32, 33 ; however, they display significant differences in single channel conductance, phar- 
macological profiles 32, 33 , and kinase activity 17 . Unlike TRPM7 whose channel activity is known to be controlled by 
PIP 2 34 , how TRPM6 is gated is not clear. 

Here, we show that TRPM6 channel activity and TRPM6-mediated Mg 2+ influx are controlled by PIP 2 levels. 
Depletion of PIP 2 by G q -linked receptor activation, by depolarization-induced activation of voltage-dependent 
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phosphatase (Ci-VSP), and by chemical translocation of 5-phospha- 
tase, can all efficiently inactivate TRPM6. Neutralization of the posi- 
tively charged residues in the TRP domain leads to dysfunctional or 
nonfunctional mutants with reduced single channel activity by PIP 2 , 
suggesting that these positively charged residues are likely to be the 
putative PIP 2 binding sites. Furthermore, we demonstrate that 
the kinase domain of TRPM 6 interacts with PLC isoforms, although 
the interaction is not necessary for PLC-induced TRPM6 channel 
inactivation. These results indicate that PIP 2 controls TRPM6 gating, 
and perhaps Mg 2+ homeostasis, under various physiological/patho- 
logical conditions. 

Results 

Phospholipase C stimulation inactivates the TRPM6 channel. 

Magnesium homoestasis is tightly controlled by many hormone/ 
receptor interactions. In order to understand how TRPM6 gating 
is regulated, we first determined whether activation of G protein 
coupled receptor induced PIP 2 hydrolysis has any effect on 
TRPM6 activation. TRPM6 was transfected to HEK-293 cells stably 
expressing the Ml receptor (HM1). Under whole cell configuration, 
TRPM6 current was small right after rupture, and ran up with time 
when intracellular free Mg 2+ concentration was decreased after 
pipette solution dialyzed into the cell. Upon TRPM6 current reach- 
ing a steady-state, 200 uM carbachol (CCh) was applied to the cell 
(Fig. 1A). TRPM6 current was rapidly and almost completely inhi- 
bited by CCh application (Fig. 1A-B), suggesting that activation 
of the Ml receptor by CCh inactivates TRPM6 channel activity. In 
agreement with this notion, CCh failed to inhibit TRPM6 current in 
HEK-293 cells transfected with TRPM6 but without over-expression 
of Ml receptor (sFig. 1). Ca 2+ release induced by CCh was detected in 
HM1 cells but not in HEK-293 cells, further suggesting that HEK- 
293 cells lack Ml receptor (sFig. 1). In order to determine whether 
Ml receptor stimulation induced inactivation of TRPM6 is through 
the PLC pathway, we tested the effect of PLC inhibitor U-73122. 
However, we found that both the PLC inhibitor U-73122 and its 
inactive isoform U-73343 inhibited TRPM6 and TRPM7 currents 
(sFig. 2). The inhibitory effects of U-73122 could be caused by direct 
inhibition of the channels, or through the partial agonist effect of 
U-73122 on PLC 35 . 

Although HM1 cells have much lower expression of endogenous 
TRPM7 (~ 200 pA) 36 , to ensure that the CCh inhibited current is 
indeed TRPM6 current, we applied 200 uM 2-APB before applica- 
tion of CCh. We have previously demonstrated that 200 uM 2-APB 
can completely inhibit TRPM7 currents, but fully potentiate TRPM6 
currents 33 . As shown in Fig. 1C, 2-APB at 200 uM markedly 
enhanced current amplitude by 100%, indicating that the currents 
recorded in HM1 cells that were largely inhibited by CCh were 
indeed TRPM6 currents. 

Stimulation of the Ml receptor by CCh can activate PLC-f3 which 
hydrolyzes PIP 2 and generates diacylglycerol (DAG) and inositol 
1,4,5-trisphosphate (IP 3 ). We hypothesized that depletion of PIP 2 
inactivates TRPM6 channel activity, similar to the effects of PIP 2 
on TRPM7 34 . Indeed, inclusion of DAG or IP 3 in the pipette solution 
did not influence TRPM6 current amplitude (data not shown). 
Furthermore, we applied GTPyS in the pipette solution to mimic 
the effect of Ml receptor activation by CCh. Consistent with the 
effects of CCh on TRPM6 current in HM1 cells, GTPyS inhibited 
TRPM6 currents in a dose-dependent manner (Figure 1E-F). These 
results indicate that Gq linked receptor activation can inhibit 
TRPM6 channel activity. Since GTPyS activates several downstream 
pathways, we further tested the effects of activators of PKC and 
cAMP on TRPM6 currents. As shown in sFig. 3, application of 
5 uM Foskolin to activate cAMP did not change TRPM6 current 
amplitude, consistent with previous report 37 . Perfusing the cells with 
PDBu to activate PKC could not produce significant changes on 
TRPM6 currents either. Furthermore, GTPyS included in the pipette 
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Figure 1 | TRPM6 current is inhibited by CCh stimulation of the Ml 
receptor. (A) A representative recording in HM 1 cells transfected with 
TRPM6. Currents were elicited by 100 ms voltage ramps ranging from 
-100 to +100 mV. Application of 200 uM CCh to activate Ml receptor 
completely inhibited TRPM6 current. (B) Time-dependent changes of 
outward TRPM6 current measured at +100 mV before and after 200 uM 
CCh application. (C) Potentiation of the currents by 2-APB was used to 
confirm that the recordings were TRPM6 currents. (D) Mean current 
densities of TRPM6 before and after CCh (n= 10, * p<0.05). (E) 
Representative traces of TRPM6 currents recorded from different cells 
using pipette solution with or without 1 mM GTPyS. (F) Concentration- 
dependent effects of GTPyS on TRPM6. Currents recorded at various 
GTPyS concentrations were normalized to the current recorded without 
GTPyS in the pipette solution. Best fit of the dose-response curve yielded 
IC 5O =10.1 uM (n=5-9 cells for each GTPyS concentration). 

solution could indeed cause PIP 2 deletion, as detected by the trans- 
location of the PH-GFP, an indicator of PIP 2 . (sFig. 3C-D). These 
results suggest that GTPyS elicited inactivation of TRPM6 is mainly 
through PLC activation induced PIP 2 hydrolysis. 

Direct removal of PIP 2 by chemical translocation of 5-phosphatase 
for PtdIns(4,5)P 2 to the plasma membrane inactivates TRPM6. To 

determine whether depletion of PIP 2 underlies the mechanism by 
which TRPM6 is inhibited by Ml receptor activation, we first applied 
a chemical dimerization system induced by Rapamycin (Rap) to 
deplete PIP 2 . The membrane-localized Rap binding protein PM- 
FRB and the fluorescent cytoplasmic enzyme construct FKBP-5-phos- 
phatase (5-ptase) were co- transfected with TRPM6 into HM1 cells. 
Application of the chemical inducer Rap will recruit 5-ptase to the 
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Figure 2 | Membrane tethering of type IV 5-ptase reduces membrane PIP 2 and suppresses TRPM6 current. (A) Schematic diagram demonstrating 
how the rapamycin (Rap) inducible PIP 2 specific phosphatase system works. (B) Representative recordings of TRPM6 before and after Rap application 
in cells transfected with GFP-TRPM6, FKBP-Inp54p and Lyn-FRB. (C) Time course of TRPM6 inhibition by perfusing 10 jiM Rap. (D) Average 
inhibition of TRPM6 by 10 uM Rap (n= 12) in cells transfected with TRPM6, FKBP-Inp54p and Lyn-FRB. (E) Representative recordings of TRPM6 
before and after Rap application in cells transfected with TRPM6 alone. (F) Time-dependent changes of TRPM6 before and after 10 uM Rap in cells shown 
in (E). (G) Mean current densities of TRPM6 before and after Rap. Rap did not produce inhibition on TRPM6 in cells transfected with TRPM6 alone. 



plasma membrane, thereby rapidly and irreversibly convert PIP 2 to 
PI(4)P by removing the 5' position phosphate at the triphosphoino- 
sitol ring (Fig. 2A) 38, 39 . As shown in Fig. 2, in the cells co-expressing 
TRPM6 with 5-ptase and the anchor protein FBR, application of 
10 pM Rap substantially reduced TRPM6 current amplitude 
(Fig. 2B-C) with an average inhibition of 90%± 9% (Fig. 2D). Rap 
(10 (J.M) did not produce noticeable effects on TRPM6 current in the 
cells without over-expression of 5-ptase and anchor protein FRB 
(Fig. 2E-G). Moreover, when TRPM6 was co-transfected with FRB 
and the mutant 5-ptase (D281A) 38 , application of 10 pM Rap failed to 
reduce TRPM6 current amplitude (sFig. 4). Thus, the results in Fig. 2 
demonstrate that activation of 5-ptase to deplete PIP 2 can fully 
inactivate TRPM6 channels, providing strong evidence that PIP 2 but 
not PI(4)P is required for TRPM6 channel function. 

Simultaneous changes of PIP 2 level and TRPM6 channel activity. 

The membrane PIP 2 levels can be modulated by a variety of stimuli 
and intracellular events. A recently identified voltage-sensitive phos- 
phatase from Ciona intestinalis (Ci-VSP) possesses voltage-sensitive 
phosphatase activity 40 . Ci-VSP can convert PIP 2 to PI(4)P upon 
depolarization, but does not have phosphatase activities at hyper- 
polarized voltage 40 ' 41 . Co-transfection of TRPM6 and Ci-VSP 



provides us an opportunity to change membrane PIP 2 levels using 
TRPM6 recording protocol (voltage ramp) without applying addi- 
tional chemicals. Cells were held at —60 mV for 3 min to allow 
the pipette solution to fully dialyze into the cells. Then, a standard 
voltage ramp protocol ranging from —120 mV to +100 mV was 
used to activate Ci-VSP and to record TRPM6 currents (Fig. 3A). 
As shown in Fig. 3, TRPM6 current was rapidly decreased in the cells 
co-expressing TRPM6 and Ci-VSP (Fig. 3B, D-E); whereas, in 
cells transfected with TRPM6 alone, TRPM6 current remained 
unchanged (Fig. 3C-D, F), indicating that depolarization induced 
dephosphorylation of PIP 2 by Ci-VSP inhibited TRPM6 channel 
activity. These results provide further evidence that PIP 2 hydrolysis 
inactivates TRPM6. 

To determine the correlation of PIP 2 hydrolysis and TRPM6 chan- 
nel activity, we simultaneously monitored the changes of membrane 
PIP 2 levels and TRPM6 current amplitude. A PtdIns(4,5)P 2 reporter, 
PLCSl-PH-GFP construct, was co-transfected with TRPM6 and 
Ci-VSP to the cells. PH-GFP is localized in the cell membrane by 
binding to PIP 2 , and can be released to the cytosol when PIP 2 is 
hydrolyzed. We monitored the translocation of PH-GFP and 
recorded TRPM6 current simultaneously in the same cells. In the 
co-transfected cells, TRPM6 current was rapidly inhibited in the first 



SCIENTIFIC REPORTS | 1 : 146 | DOI: 1 0. 1 038/srep00146 



3 



www 



.nature.com/ scientificreports 




Figure 3 | Simultaneous monitoring of PIP 2 depletion and TRPM6 
inactivation. (A) The protocol for activation of Ci-VSP and recording of 
TRPM6 currents. Cells were held at — 60 mV to allow dialysis of pipette 
solution without activation of Ci-VSP. Voltage Ramp ranging from — 120 
to + 100 mV was applied to activate Ci-VSP and record TRPM6 currents. 
(B-C) Representative traces recorded before and after activation of Ci-VSP 
in cells transfected with TRPM6 and Ci-VSP (B) and TRPM6 alone (C). 
(D) Time-dependent changes of TRPM6 outward current measured at 
+ 100 mV in cells transfected with TRPM6+Ci-VSP and TRPM6 alone. 
(E-F) Normalized currents before and after depolarization from cells with 
TRPM6+Ci-VSP transfection (E) and with TRPM6 transfection alone (F) 
(n=8 for each group). (G) Representative GFP-PH domain sub-cellular 
location at different time points upon depolarization. Cells were 
co-transfected with GFP-PH, Ci-VSP and TRPM6. (H) Kinetics of cytosol 
fluorescence increment and TRPM6 current inhibition. The intensity of 
fluorescence and current amplitude of TRPM6 at each time-point were 
normalized to the initial values, respectively. The time constants obtained 
by mono-exponential fit were 24.2 ± 2.5 s for PH-GFP translocation, and 
36.5±6.9 s for TRPM6 inactivation. 

20 s, and almost completely inhibited at 200 s; whereas, in the 
TRPM6 only transfected cells, current amplitude was barely changed 
over the time course of 200 s recording (Fig. 3D). Meanwhile, the 
membrane fluorescence representing translocation of PH-GFP sig- 
nificantly decreased in the first 20 s, and slowly disappeared over the 
time course of 200 s (Fig. 3G-H). The time-dependent inhibition of 
TRPM6 current amplitude exhibited a similar kinetics to that of 



translocation of PH-GFP. The time constant was 36.5 ±6.9 s for 
TRPM6 inhibition, and for PH-GFP translocation was 24.2 ± 2.5 s 
(Fig. 3H). The kinetics of PIP 2 hydrolysis by activation of Ci-VSP 
was slower than previously reported 42 , probably due to the ramp 
protocols used in our experiments. Nevertheless, the well-correlated 
kinetics of TRPM6 inactivation and PIP 2 hydrolysis (Fig. 3H) indi- 
cates that TRPM6 channel inactivation is caused by the decrease of 
PIP 2 levels. Indeed, including water soluble PIP 2 (DiC8-PIP 2 ) in the 
pipette solution significantly reduced the inactivation speed of 
TRPM6 caused by PIP 2 hydrolysis in the TRPM6 and Ci-VSP co- 
transfected cells (sFig. 5), supporting the notion that PIP 2 is neces- 
sary for TRPM6 channel activation. 

Depletion of PIP 2 inactivates TRPM6 channel activity in perfo- 
rated-patch. Since the above results were obtained by whole-cell 
configuration, to eliminate any interference of pipette solution to 
the effects of PIP 2 on TRPM6 activity, we applied perforated-patch 
to determine the effects of PIP 2 hydrolysis on TRPM6. As shown in 
sFig. 6, in the cells co-transfected with TRPM6 and Ci-VSP, there was 
substantial TRPM6 current recorded by perforated-patch, albeit the 
current amplitude being smaller than that of whole-cell current 
amplitude. TRPM6 current was readily reduced by depolarization 
which activated Ci-VSP to hydrolyze PIP 2 , further confirming that 
PIP 2 is essential for TRPM6 channel activity. The effects of deletion 
of PIP 2 in perforated patches were also observed in TRPM7 trans- 
fected cells (sFig. 7). 

Putative PIP 2 interacting residues of TRPM6. PIP 2 regulates ion 
channel function by binding to specific residues 43 . To understand the 
mechanism by which PIP 2 activates TRPM6, we mutated the posi- 
tively charged residues equivalent to those of TRPM8, which have 
been shown to be responsible for PIP 2 modulation 44 . We created 
single and double mutants as shown in Fig. 4, and tested channel 
function of these mutants. The current amplitude of K1085Q and 
R1088Q mutants was much smaller than that of WT TRPM6 
(Fig. 4B, E), yet they can be potentiated by 2-APB (Fig. 4C-D), 
indicating that the recorded currents from K1085Q and R1088Q 
transfected cells were TRPM6. The reduced current amplitude of 
the mutants could be the consequence of disrupted interaction of 
PIP 2 with the putative binding residues. Interestingly, the mutation 
of Lys to Glu (K1098Q) completely abolished TRPM6 channel acti- 
vity (Fig. 4E), suggesting that Lys 1098 may be more crucial for TRPM6 
channel activity than Lys 1085 and Arg 1088 . The loss-of-function 
mutant K1098Q has no response to 2-APB stimulation (Fig. 4F). 
Furthermore, although both K1085Q and R1088Q mutants pro- 
duced small currents, double mutation by neutralizing both Lys 1085 
and Arg 1088 abolished channel function completely (Fig. 4E-F). Like 
K1098Q, the non-functional double mutant K1085Q/R1088Q could 
not be reactivated or potentiated by 2-APB, suggesting that the 
effects of 2-APB on TRPM6 depend on the presence of PIP 2 and 
the availability of the interaction between PIP 2 and TRPM6. This 
hypothesis is supported by the results that 2-APB could not poten- 
tiate TRPM6 after the channel is completely inactivated by PIP 2 
depletion induced by Ci-VSP activation (sFig. 8). 

To determine whether mutations disrupt plasma membrane 
localization leading to dysfunctional channels, we measured plasma 
membrane protein and total protein levels of TRPM6 and the 
mutants (Fig. 5). GADPH was used as a negative control for the 
plasma membrane protein. The ratios of plasma membrane/total 
protein of K1085Q, R1088Q, K1098Q, and K1085Q/R1088Q were 
not significantly different from that of WT TRPM6 (Fig. 5B), indi- 
cating that mutations of the positively charged residues in the TRP 
domain, but not the expression or trafficking of the channel proteins, 
disrupted TRPM6 channel function. 

PIP 2 regulates single channel activities in excised patches of 
TRPM6 and its mutant. To examine directly the effects of PIP 2 on 
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Figure 4 | PIP 2 binding sites of TRPM6. (A) Alignment of the TRP domain of TRPM6, TRPM7 and TRPM8. The highlighted residues in TRPM8 
are the PIP 2 binding sites. (B) Representative traces of WT-TRPM6 and mutants K1085Q, R1088Q, and K1098Q. (C-D) Effects of 2-APB on the mutants 
K1085Q and R1088Q. (E) Normalized mean current density of TRPM6 mutants in comparison with WT TRPM6 (n= 10 *p<0.05; ** P<0.01). Mutants 
K1098Q and the double mutants K1085Q/R1088Q did not produce any current. (F) Average changes of current amplitude by 2-APB (200 uM). 2-APB 
did not have any effect on the non-functional mutants K1098Q and K1085Q/R1088Q. 



channel activity, TRPM6 single channel currents were recorded at 
— 50 mV in inside-out patches in the DVF solution as shown in 
Fig. 6A. Similar single channel currents were not observed in 
mock-transfected cells. The slope conductance calculated based 
on single channel current amplitude at —50 mV (3.99±0.3pA), 
-80 mV (6.7±0.8pA), and -120 mV (9.8±l.lpA) was 82.7 pS 
(n=4), which is indistinguishable from the signature conductance 
of TRPM6 (83.6 pS) 33 . To further confirm that the recorded currents 
were TRPM6 currents, we applied 2-APB to the excised patches. 
Similar to the effects of 2-APB on macroscopic TRPM6 currents, 
we observed that 2-APB increased TRPM6 channel activity (Data 
not shown). Like other PIP 2 activated channels such as TRPM7 34 and 
TRPM8 44 , single channel activities of TRPM6 rundown in excised 
patches (Fig. 6A), presumably caused by depletion of PIP 2 34,44 . After 
rundown, application of PIP 2 (DiC8-PIP 2 ) recovered TRPM6 chan- 
nel activity by almost 80% (Fig. 6A, D), indicating that PIP 2 is suf- 
ficient to activate TRPM6. In contrast to the WT TRPM6 single 
channel activity, the putative PIP 2 binding residue mutant R1088Q 
exhibited extremely low open probability (Fig. 6B-C), suggesting 



that neutralization of the positive charges of the putative PIP 2 
binding sites might have impaired channel activity. Moreover, after 
rundown, application of 10 or 50 uM PIP 2 to the excised patches of 
R1088Q only recovered the channel activity by 10%, suggesting that 
the activation of R1088Q by PIP 2 is significantly decreased in com- 
parison with WT TRPM6. These results imply that R1088 is likely to 
be involved in the interaction of TRPM6 and PIP 2 . 

The putative PIP 2 binding residues are essential TRPM7 channel 
activity. TRPM6 and TRPM7 are the closest homologues with a num- 
ber of similar and distinctive properties 17, 18, 33, 45 . Previous studies have 
shown that depletion of PIP 2 inactivates TRPM7, yet the PIP 2 binding 
sites were not determined 34, 46 . Thus, we mutated the positively charged 
residues equivalent to those of TRPM6 in the TRP domain (Fig. 4A). 
Unlike TRPM6, all the single mutants of TRPM7 exhibited channel 
activity (Fig. 7A-B), although the current amplitude of K1112Q was 
significantly smaller than that of WT TRPM7. The double-mutant 
K1112Q/R1115Q and the triple mutant K1112Q/R1115Q/K1125Q, 
however, completely lost channel function. Moreover, current 
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amplitude of the cells transfected with K1112Q/R1115Q or K1112Q/ 
Rl 1 1 5Q/K1 125Q (Fig. 7A) was even smaller than that of mock-trans- 
fected cells, suggesting that the loss-of-function mutants K1112Q/ 
R1115Q and K1112Q/R1115Q/K1125Q behave as a dominant- 
negative to the endogenous TRPM7 channel. Given that PIP 2 has 
been shown to be necessary and sufficient for TRPM7 channel 
activation, it is likely that the conserved positively charged residues 
in the proximity of the TRP domain are essential for PIP 2 and 
TRPM7 interaction. 

Since PIP 2 is the only activator for TRPM7, to further confirm that 
double mutant and the triple mutant of TRPM7 are indeed loss-of- 
function mutants, we tested the effects of NH 4 C1 on these mutants. 
Perfusion of NH 4 C1 has been proposed to increase intracellular 
pH thereby releasing PIP 2 sequestered by other cations such as H + , 
Mg 2+ , or other polyvalent cations 46 . For these experiments, we 
included 1 mM free Mg 2+ in the pipette solution. As shown in 
Fig. 7C, after perfusion with 30 mM NH 4 C1, Kl 1 12Q current ampli- 
tude was increased by 5 -fold; whereas the triple mutant had no 
response to NH 4 C1, suggesting that the double mutant (data not 
shown) and the triple mutant are nonfunctional mutants. 

PIP 2 depletion abolishes Mg 2+ currents and Mg 2+ influx through 
TRPM6. Next we examined whether PIP 2 depletion affects Mg 2+ 
influx through TRPM6. As shown in Fig. 8 A, Mg 2+ currents recorded 
using isotonic Mg 2+ external solution were significantly increased 
by 200 uM 2-APB. After depletion of PIP 2 in cells co-transfected 
with TRPM6 and Ci-VSP, the Mg 2+ current was fully inhibited 
(Fig. 8A-B). Similar results were observed using CCh and GTPys 



to delete PIP 2 (data not shown). We further measured Mg 2+ influx 
through TRPM6 before and after PIP 2 depletion. The membrane 
permeable Mg 2+ indicator KMG104-AM was used to measure 
changes of intracellular Mg 2+ induced by perfusing cells with iso- 
tonic Mg 2+ solution. As shown in Fig. 8 C, significant increase in the 
normalized fluorescence intensity (F/F 0 ) was induced by perfusion of 
isotonic Mg 2+ in TRPM6-transfected cells in comparison with the 
non-transfected cells. Moreover, the change of Mg 2+ influx was 
largely diminished in the cells perfused with 200 (xM CCh in DVF 
and isotonic Mg 2+ solutions to hydrolyze PIP 2 . These results indicate 
that hydrolysis of PIP 2 can induce TRPM6 inactivation and abolish 
TRPM6 mediated Mg 2+ influx, suggesting that activation of PLC by 
Gq-linked receptors may play an pivotal role in regulating Mg 2+ 
homeostasis by controlling TRPM6 channel activity. 

TRPM6 kinase domain interacts with PLC, but the interaction is 
not required for PLC induced TRPM6 inactivation. As shown in 
previous studies, the TRPM7 channel kinase domain interacts with the 
C2 domain of PLC, especially PLC- [32 isoform 34 . Since TRPM6 and 
TRPM7 share high homology, we tested whether the TRPM6 kinase 
domain binds to the PLC isoforms. As shown in sFig. 9, the TRPM6 
kinase domain primarily binds to the PLC P isoforms, including 
PLCpi-3. Unlike GST-M7 KIN, GST-M6 KIN did not pull-down 
the PLC-Sl. This result suggests that the kinase domain of TRPM6 
and TRPM7 preferentially interact with various PLC isoforms. 

To understand whether interaction of the kinase domain with PLC 
is necessary for PLC activation induced inhibition of TRPM6 and 
TRPM7, we made truncation mutants of TRPM6 (M6-Akinase) 
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Figure 6 | Effects of PIP 2 on channel activities of TRPM6 and its mutant in the inside-out patches. (A) Representative recordings of WT TRPM6 single 
channel currents in an inside-out patch at -50 mV recorded in divalent free solution (DVF). PIP2 (DiC8-PIP 2 ) at 10 uM was applied to the patch after 
completely rundown of channel opening. (B) Single channel currents of R1088Q recorded at —50 mV right after formation of excised patch, after 
rundown, and upon application of 10 uM PIP 2 . (C) Average open probability of TRPM6 (n=8) andR1088Q (n=7) before and after completely rundown. 
(D) Percentage of channel open probabilities rescued by 10 uM PIP 2 after rundown of TRPM6 and R1088Q. The effect of DiC8-PIP 2 was reversible and 
reproducible in the same patches. 



and TRPM7 (M6-Akinase) by deleting the kinase domain at the 
C-terminus. Currents recorded from M6-Akinase were verified by 
response to 2-APB (sFig. 9 C). Since we have shown that GTPyS 
inactivates TRPM6 mainly through activation of PLC pathway, we 
used GTPyS to activate PLC. As shown in sFig. 9 C-F, inclusion of 
GTPyS in the pipette solution markedly inhibited M6-Akinase and 
M7-Akinase currents, indicating that the interaction between kinase 
domain and the PLC is not necessary for the PLC activation induced 
TRPM6 and TRPM7 channel inactivation. Similar results were 
also obtained using CCh to induce inactivation of M6-Akinase 
and M7-Akinase expressed in the HM1 cells. CCh (200 uM) com- 
pletely inhibited M6-Akinase and M7-Akinase currents (n=5 for 
each group). 

Discussion 

We reported several new findings here. First, we demonstrated that 
PIP 2 is necessary and sufficient for TRPM6 activation. Second, we 
identified that the positively charged residues in the TRP domain are 
likely to participate in TRPM6 and PIP 2 interactions. Third, we 
found that TRPM6 interacts with PLC isoforms, but the interaction 
is not necessary for PLC induced TRPM6 or TRPM7 inactivation. 
Fourth, inactivation of TRPM6 by PIP 2 depletion eliminates Mg 2+ 
currents and Mg 2+ influx, suggesting that PIP 2 regulates Mg 2+ home- 
ostasis by controlling TRPM6 gating. Our results indicate that 
regulation of TRPM6 by PIP 2 through PLC activation pathway 
may serve as an important pathway to gate TRPM6 and control 
Mg 2+ homeostasis. 

TRPM6 plays a vital role in Mg 2+ homeostasis, as demonstrated by 
the mutations of TRPM6 that cause HSH 5 - 10 . Although it has been 
reported that TRPM6 expression can be regulated by various factors, 
including estrogen 26- 47 , Angll 27 , acidosis/alkalosis conditions 28 , 
immunosuppressant tacrolimus 29 , diuretics Thiazide 30 , and EGF 31 , 



how TRPM6 is gated has remained elusive. Here we provided several 
lines of evidence demonstrating that TRPM6 gating is controlled by 
cellular PIP 2 levels. Activation of PLC-coupled Ml receptor by CCh 
completely suppressed TRPM6 current amplitude (Fig. 1). This 
result is supported by pharmacological activation of PLC induced 
suppression of TRPM6 (Fig. 1E-F). Using the recently developed 
approach for the study of phosphoinositide signaling, which is based 
on the rapamycin-induced heterodimerization of the rapamycin 
(Rap) binding domain of mTOR (FRB) and FKBP12 (Fig. 2) 38 ' 3S , 
we showed that application of Rap rapidly and fully inactivated 
TRPM6 channel activity (Fig. 2). Rap induced translocation of 
over-expressed type-IV 5-phosphatase (5-Ptase) specifically hydro- 
lyzes PIP 2 without involvement of other signaling factors. Thus, Rap 
induced inactivation of TRPM6 strongly indicates that PIP 2 is 
required for TRPM6 channel function. This conclusion was further 
supported by the fact that TRPM6 can be inactivated by Ci-VSP 
elicited PIP 2 depletion (Fig. 3). Furthermore, application of PIP 2 
in excised patches can reactivate TRPM6 after rundown (Fig. 6). 
Thus, it appears that PIP 2 is necessary and sufficiently for TRPM6 
channel activity. 

PIP 2 regulates a large number of ion channel functions by inter- 
acting with specific residues or domains 43 . However, due to the lack 
of crystal structure information of most of ion channels, putative 
PIP 2 interacting sites are usually predicted based on direct binding 
experiments with PIP 2 , or changes of PIP 2 affinity evaluated by 
channel functions between Wide-Type and mutants 43 . Positively 
charged residues at both N- and C-termini have been shown to be 
putative PIP 2 interacting domains 43, 48 ~ 50 . Rohacs and colleagues 
demonstrate that the basic residues within the TRP domain are likely 
the putative PIP 2 binding sites for TRPM8, TRPV5 and TRPM5 44 . 
Nilius and colleagues found that a PH-like domain at the more distal 
C-terminal region of TRPM4 may be involved in PIP 2 -dependent 
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Figure 7 | Mutation of the putative PIP 2 binding residues of TRPM7 
disrupts channel function. (A) Representative traces of WT TRPM7 and 
its mutants. The double mutant K1112Q/R1115Q (DM) and the triple 
mutant K1112Q/R1115Q/K1125Q did not produce currents. (B) Mean 
current densities of the WT and its mutants (n=8-13, *p<0.05). (C) 
Effects of NH 4 C1 on Kl 1 12Q, and triple mutant Kl 1 12Q/R1 1 15Q/K1 125Q 
(TM) were tested by using pipette solution containing 1 mM free Mg 2+ . 
Perfusion of 30 mM NH 4 C1 can increase intracellular pH thereby releasing 
PIP 2 sequestered by Mg 2+ , thus, the effects of PIP 2 on the channel activity is 
manifested. NH 4 C1 significantly increase the current amplitude of 
Kl 1 12Q, but failed to induce any change in the DM and TM transfected 
cells, indicating that the double mutant K1112Q/R1115Q and the triple 
mutant Kl 1 12Q/R1 1 15Q/K1 125Q have no ability to sense the changes of 
PIP, levels. 



regulation 51 . We mutated the conserved basic residues at the TRP 
domain of TRPM6, which are equivalent to the putative PIP 2 -bind- 
ing sites of TRPM8 44 . The mutants K1085Q and R1088Q displayed 
largely diminished channel activity, whereas K1098Q and the double 
mutant K1085Q/R1088Q completely lost channel function. Since 
PIP 2 is the only activator of TRPM6, the decreased channel function 
of the mutants suggests that the mutations disrupt TRPM6 activation 
by PIP 2 . Moreover, it seems that mutant R1088Q has reduced activity 
by PIP 2 , as evidenced that PIP 2 at 10 (J.M recovered channel 
activity by 80% of WT, but only 10% of R1088 (Fig. 6). Mutation 
of the conserved residues of TRPM7 also markedly reduced channel 



activity. Therefore, it is likely that the basic residues at the TRP 
domain may participate in channel-PIP 2 interactions. Further 
studies are required to confirm whether the positively charged resi- 
dues at the TRP domain are PIP 2 binding sites. Excitingly, Hansen 
and colleagues have recently discovered the structure basis of inter- 
action between PIP 2 and Kir2.2 52 . Putative PIP 2 binding residues at 
the cytoplasmic domain predicted by previous mutagenesis studies 43 
are proven to interact with PIP 2 . A conserved non-specific phospho- 
lipid binding region is also shown to constitute the PIP 2 -binding 
interface with the cyloplasmic domain for PIP 2 binding. Thus, it is 
conceivable that the predicted PIP 2 binding residues at the TRP 
domain may be proven to be the PIP 2 binding sites in the future. 

Like TRPM7 34 , TRPM6 kinase domain interacts with different 
PLC isoforms (sFig. 9). However, we found that the interaction of 
M6-KIN and M7-KIN with PLCs is not required for modulation of 
TRPM6 and TRPM7 by PLC activation, as evidenced by the fact that 
the truncation mutants of TRPM6 and TRPM7 without the kinase 
domain can also be inactivated by PLC activation (sFig. 9). Although 
it is not necessary for PLC-induced regulation on TRPM6 and 
TRPM7, given the different distribution profiles of the PLC iso- 
forms 53- 54 , the interaction or co-localization of TRPM6 and 
TRPM7 with different isoforms of PLC may underlie different, yet 
to be determined physiological/pathological functions. 

Our results seem to indicate that TRPM6 is exclusively gated by 
PIP 2 . The loss-of-function mutants of putative PIP 2 binding sites 
cannot be potentiated by 2-APB (Fig 4). Similarly, inactivated 
TRPM6 by completely deletion of PIP 2 has no response to 2-APB 
(sFig. 8). PIP 2 may also be the exclusive activator of TRPM7 34 . The 
important role of PIP 2 on TRPM7 activation was demonstrated in 
2002 34 , and the conclusion was supported by other studies 46, 55 . 
However, Takezawa and colleagues reported that activation of PLC 
through endogenous Ml-muscarinic receptors has little effect on 
TRPM7 over-expressed in HEK cells 56 . Rather, they suggested that 
cAMP enhances TRPM7 activity. It was reported later that failure of 
inhibition of TRPM7 by CCh using endogenous Ml receptor 56 was 
probably due to the minor PIP 2 breakdown by limited endogenous 
Ml receptor 57 . The effects of cAMP on TRPM7 activation 56 was not 
observed by Langeslag and colleagues 57 . Indeed, Langeslag and 
colleagues elegantly demonstrated a positive correlation between 
TRPM7 channel activity and PIP 2 level 57 , in agreement with the 
results that TRPM7 activity requires PIP 2 as shown in Fig. 7 and as 
previous reported 34 . However, under perforated-patch configura- 
tion, Langeslag and colleagues observed that Bradykinin induced a 
transient inhibition of TRPM7 currents 57 , because PIP 2 levels were 
quickly recovered in N1E-115 cells. They therefore suggested that 
PLC-coupled agonists activates TRPM7 57 . Although it was unclear 
why PIP 2 levels were quickly recovered under the perforated-patch 
configuration 57 , it is plausible that other signaling pathways were 
involved in PIP 2 synthesis under their experimental conditions. 
Nevertheless, using perforated patch configuration, we demon- 
strated that activation of the voltage-sensitive Ci-VSP by depolariza- 
tion can completely inhibit both TRPM6 (sFig. 6) and TRPM7 
(sFig. 7) channel activities. Therefore, it seems that PIP 2 is a strong 
and exclusive activator of TRPM6 and TRPM7. 

Mg 2+ permeation through TRPM6 was substantially enhanced 
when TRPM6 channel activity was increased by 2-APB. However, 
after depletion of PIP 2 by CCh, Mg 2+ current and Mg 2+ influx 
were largely blocked (Fig. 8). This PLC-activation induced Mg 2+ 
regulation may produce physiological/pathological significance. 
In the distal convoluted tubule (DCT) where TRPM6 is highly 
expressed, Ca 2+ /Mg 2+ -sensing receptors (CaSR) are also abundantly 
expressed 58 . Stimulation of CaSR activates PLC, especially the (3 and y 
isoforms 59 , which therefore may influence the magnesium re-absorp- 
tion by regulation of TRPM6 channel function. Several peptide 
hormones such as PTH, calcitonin, glucagon, and AVP enhance 
magnesium reabsorption in the DCT 60 . Groenestege and colleagues 
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Figure 8 | Depletion of PIP 2 inhibits Mg 2+ currents and eliminates Mg 2+ 
influx through TRPM6. (A) Mg 2+ currents recorded using isotonic 
extracellular Mg 2+ in cells transfected with either TRPM6 alone (blue) or 
TRPM6+Q-VSP (red). Mg 2+ current (blue) was significantly increased by 
2-APB (green), and completely inhibited after PIP 2 depletion induced 
by activation of Ci-VSP. (B) Normalized Mg 2+ current amplitude (n=9). 

(C) Mg 2+ influx through TRPM6 under control conditions and after PIP 2 
hydrolysis by 200 uM CCh. Changes of Mg 2+ influx were measured in 
isotonic Mg 2+ extracellular solution. Note the normalized fluorescence 
intensity (F/F 0 ) induced by isotonic Mg 2+ was largely diminished by CCH. 

(D) Average changes of F/F 0 (n=33~36 in each group). 

very nicely showed that EGF regulates Mg 2+ re-absorption by regu- 
lation of TRPM6 trafficking to the cell membrane 37 ' 47 . A mutation of 
the EGF gene encoding pro-EGF causes renal Mg 2 " 1 " wasting 37, 47 . As 
EGF can also activate PLC-y, leading to PIP 2 depletion, there might a 
fine-tuned regulation regarding how EGF modulates TRPM6. 
Nonetheless, it will be of great interest to investigate whether and 
how the other Mg 2+ regulating hormones may regulate the Mg 2+ 
homeostasis through Mg 2+ gate-keeper TRPM6. 

In conclusion, we have demonstrated at various levels that PIP 2 is 
required for TRPM6 channel activity, and that disruption of the 
putative PIP 2 binding sites results in dysfunctional or non-functional 
TRPM6 channels. As TRPM6 is considered as a Mg 2+ gate-keeper, 
understanding the critical role of PIP 2 in TRPM6 channel function 
may facilitate our understanding of how Mg 2+ contributes to physio- 
logical/pathological processes, as well as the development of thera- 
peutic approaches for hypomagnesemia treatment. 

Methods 

Molecular biology. TRPM6 plasmid (pCINeo/IRES-GFP) was kindly provided by 
Dr. Joost G. J. Hoenderop. Mutations of TRPM6 and TRPM7 were made by using the 
QuickChange Site-Directed Mutagenesis Kit (Stratagene) following the manufac- 
ture's instruction. The primers sequences will be available upon request. 

Electrophysiology. Whole-cell and single- channel currents were recorded using an 
Axopatch 200B amplifier. Data were digitized at 5 or 10 kHz, and digitally filtered off- 
line at 1 kHz. Patch electrodes were pulled from borosilicate glass and fire-polished to 
a resistance of ~-3 Mfl when filled with internal solutions. Series resistance (Rs) was 
compensated up to 90% to reduce series resistance errors to <5 mV. Cells with Rs 
bigger than 10 MQ were discarded 21 . 

For whole cell current recordings, voltage stimuli lasting 250 ms were delivered at 
1 s intervals, with either voltage ramps or voltage steps ranging from — 120 to 
+ 100 mV. A fast perfusion system was used to exchange extracellular solutions, with 
complete solution exchange achieved in about 1 to 3 s 32 . The internal pipette solution 
for TRPM6 or TRPM7 whole cell current recordings contained (in mM) 145 Cs- 
methanesulfonate (CsS0 3 CH 3 ), 8 NaCl, 10 EGTA, and 10 HEPES, with pH adjusted 
to 7.2 with CsOH. Ca 2+ was adjusted to various concentrations based on calculation 



using MaxChelator 32 . The standard extracellular Tyrode's solution for whole cell 
recording contained (mM): 145 NaCl, 5 KC1, 2 CaCl 2 , 10 HEPES and 10 glucose; 
pH was adjusted to 7.4 with NaOH. Divalent-free solution (DVF) contained 
(mM) 145 NaCl, 20 HEPES, 5 EGTA, 2 EDTA and 10 glucose, with estimated free 
[Ca 2+ ]< 1 nM and free [Mg 2+ ] ~10 nM at pH 7.4. MaxChelator was used to cal- 
culate free Ca 2+ and free Mg 2+ concentrations. 

Single channel currents were recorded in inside-out patches. The extracellular DVF 
solution contained (mM) 145 Na-methanesulfonate (NaS0 3 CH 3 ), 20 HEPES, 5 
EGTA, 2 EDTA, and 10 glucose, with pH adjusted to 7.4; internal solution contained 
(mM) 145 CsS0 3 CH 3 , 8 NaCl, 10 EGTA, and 10 HEPES, with pH adjusted to 7.2. 

Rapamycin induced PIP 2 depletion. The two plasmids PM-FRB-CFP and mRFP- 
FKBP-5-ptase were kindly provided by Tamas Balla. The phosphatase dead mutant 
CF-InP-D281A was kindly provided by Dr. Thomas Mayer. PM-FRB-CFP and 
mRFP-FKBP-5-ptase were co-transfected with TRPM6 to the cells, and 10 uM 
rapamycin was applied to induce 5-ptase translocation to deplete PIP 2 . In the case of 
testing phosphatase dead mutant, PM-FRB-CFP and CF-InP-D28lA were co- 
transfected with the channel to HEK293 cells. 

Western blot. The HEK 293 cells transfected with TRPM6 were allowed 48 hours to 
express protein. Cell lysate was used to blot with anti-GFP antibody (Neuromab). The 
ECL chemiluminescence was detected and captured by FUJI LAS-3000 Imaging 
System. 

Membrane protein detection. HEK 293T cells were transfected with WT TRPM6 
(GFP-tagged) or mutants. After 48 h, cell membrane protein biotinylation was per- 
formed with Pierce Cell Surface Protein Isolation Kit (catalogue no. 89881). Briefly, 
monolayers of cells were washed with cold PBS and then incubated with Sulfo-NHS- 
SS-Biotin solution at 4 'C for 30 min. After incubation, the reaction was quenched, 
and the cells were washed 3 times with TBS solution (50 mM tris, 150 mM Nacl, 
pH 7.6), followed by lysing procedure for 30 min on ice. After centrifugation for 
10 min at 16,000 g, supernatant passed a sufficient amount of Neutr Avidin Agarose 
(Pierce Biotechnology), and the flow-through was kept as cytosol portion. Agarose 
beads with protein were washed 3 times and then boiled with IX SDS buffer con- 
taining 50 mM DTT at 95 °C, and supernatant containing released membrane protein 
was kept for analysis by Western blotting. 

Ratio Ca 2+ imaging experiments. Cells were plated on 25 mm glass coverslips 
loaded with 1 uM Fura-2 for 20 minutes. Non- incorporated dye was washed away 
using a HEPES -buffered Saline Solution (HBSS) containing (in mM) 20 HEPES, 10 
glucose, 1.2 MgCl2, 1.2 KH2P04, 4.7 KC1, 140 NaCl, 1.3 Ca 2+ (pH 7.4). Ca 2+ influx 
was measured by perfusing the cells with Tyrode's solution containing 20 mM Ca 2+ 
after the cells were perfused with nominal Ca 2+ and Mg 2+ free solution for 2 minutes. 
Ionomycin (Iono) at 1 mM was applied as internal control. Fluorescence intensities at 
510 nm with 340 nm and 380 nm excitation were collected at a rate of 1 Hz using 
CoolSNAP HQ2 (Photometries) and data were analyzed using NIS-Elements 
(Nikon) 21 . The 340 : 380 nm ratio in the presence of 20 mM Ca 2+ was normalized to 
that of 1 JJ.M Ionomycin (Iono) elicited Ca 2+ signal. 

Detection of Mg 2+ influx. HM1 cells were cultured on poly- L- lysine- coated cover 
slips in DMEM-F12 medium at 37°C. Cells were then loaded with KMG104-AM dye 
at a final concentration of 20 uM for 30 min in standard tyrode solution without 
added magnesium. Tyrode solution contained: 135 mM NaCl, 4 mM KG, 10 mM 
glucose, 10 mM HEPES, 2 mM CaCl 2 , pH 7.4. Fluorescence intensities at 510 nm 
with 488 nm excitation were collected very 10 s using CoolSNAP HQ2 
(Photometries) and data were analyzed using NIS-Elements (Nikon). F 0 refers to the 
initial fluorescence recorded at t— 0 s. Mg 2+ influx was measured by perfusing the 
cells with isotonic Mg 2+ solution containing (mM) 120 MgCl 2 , 10 Hepes, 10 glucose, 
and pH was adjusted to 7.4. 

GST-pulldown purification assay. HEK-293T cells in 10-cm dishes were transfected 
using LipofectAMINE 2000 (Life Technologies, Rockville, MD) with 1 0 ug expression 
vector containing either PLC-pl-4, PLC-yl, or GFP-PLC-51 cDNAs. After 48 h, 
cells were lysed in 2 ml of RIPA buffer (50 mM Tris-HCl at pH 7.4, 150 mM sodium 
chloride, 1 mM EDTA, 1% IGEPAL CA-630, 0.5% (w/v) deoxycholate, 0.1% (w/v) 
SDS and 10 mM iodoacetamide). GST-TRPM6 kinase, GST-TRPM7 kinase or 
GST alone was expressed and purified with glutathione-4B Sepharose beads. Bound 
GST-kinases, GST or beads were incubated with lysates containing various PLCs for 
12 h. Bound proteins were washed three times in RIPA buffer, eluted by boiling in 
2XSDS-polyacrylamide gel electrophoresis (PAGE) sample buffer and then sepa- 
rated by SDS-PAGE. PLC-pl-4 and PLC-yl were detected by western blotting with 
the following specific anti-PLC antibodies from Santa Cruz Biotechnology (Santa 
Cruz, CA): PLCpi(G-12), PLCp2(Q-15), PLCp3(C-20), PLCp4(C-18), 
PLCyl(E-12). An anti-GFP antibody (Clontech) was used to detect GFP-tagged 
PLC-51. A horseradish peroxidase -linked antibody to rabbit Ig (Amersham 
Pharmacia Biotech, Piscataway, NJ) and SuperSignal West Dura substrate was used 
for all chemiluminescent detection (Pierce, Rockford, IL). 

Data analysis. Pooled data are presented as mean ± SEM. Dose-response curves 
were fitted by an equation of the form E — £ max {l/[l + (EC 50 /C)"]}, where E is the 
effect at concentration C, E max is maximal effect, EC 50 is the concentration for half- 
maximal effect and n is the Hill coefficient 33 . Statistical comparisons were made using 
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ANALYSIS of variance (ANOVA) and two-tailed f-test with Bonferroni correction. 
P < 0.05 indicated statistical significance. 
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